LARGE DEVIATIONS FOR STATISTICS OF JACOBI PROCESS 



N. DEMNlQAND M. ZANI0 



Abstract. This paper is aimed to derive large deviations for statistics of Jacobi 
process already conjectured by M. Zani in her Thesis. To proceed, we write in a more 
simple way the Jacobi semi-group density. Being given by a bilinear sum involving 
Jacobi polynomials, it differs from Hermite and Laguerre cases by the quadratic form 
of its eigenvalues. Our attempt relies on subordinating the process using a suitable 
random time-change. This will give an analogue of Mehler formula whence we can 
recover the desired expression by inverting some Laplace transforms. Once we did, an 
adaptation of Zani's result in the non-steepness case will provide the required large 
deviations principle. 

1. Introduction 

The Jacobi process is a Markov process on [—1, 1] given by the following infinitesimal 
generator: 

for some real p, q, defined up to the first time when it hits the boundary. In fact, it 
belongs to the class of diffusions associated to some families of orthogonal polynomials, 
i.e. the infinitesmal generator admits an orthogonal polynomials basis as eigenfunctions 
([2]) such as Hermite, Laguerre and Jacobi polynomials . More precisely, if Pn denotes 
the Jacobi polynomial with parameters a, (5 > —1 defined by : 

yx,Pr-\ _ ( a + X )« n ( _ _ , _ , a , i _. , 1 . 1 



P%' p (x) = v n] ' - 2*1 ( -n, n + a + p + l,a + l; J , x e hM], 
then we can see that : 

if = - n (n + a + (3 + l)P^ 

for p = — {(3 + a + 2) and q = (3 — a. In 1964, Wong resolved the forward Kolmogorov 
or Fokker-Planck equation (see [T7], |15j ) 

dy[B(y)p] - d y [A(y)p] = d t p, p = p t {x, y), 

where B, A are polynomials of degree 2, 1 respectively, and gave the principal solution 
(po(x,y) = 8 x {y)) using the classical Sturm-Liouville theory. This gives rise to a class of 
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stationnary Markov processes satisfying : 



■'■2 



lim p t (x,y) = / W(x)p t (x,y)dx = W(y) 

where W is the density function solution of the corresponding Pearson equation (|17j). 
In our case, pt has the discrete spectral decomposition : 

(1) p t (x,y) = [Y J {Rn)- 1 e- X - t P^{x)P^{y)\ W{y), x,y G [-1,1] 

\n>0 J 

where 

X n = n(n + a + p + l), W(y) - ' ' ■~ /;) " ( ' + 



2«+/3+i J B(a + 1,(3+ 1) 
with B denoting the Beta function and (PQ, p. 99) : 

o|3_llpa,/9||2 r(a + /3 + 2) (a + l) n (/3 + l) w 

ttrU-W^n \\L*([-i,i],w(y)dy)- 2ri + a + P + ir(a + /3 + n + l)n\ 

Few years later, Gasper (|10j) showed that this bilinear sum is the transition kernel 
of a diffusion and that is a solution of the heat equation governed by a Jacobi opera- 
tor, generalizing a previous result of Bochner for ultraspherical polynomials ([6]). It is 
worthnoting that \ n has a quadratic form while in the Hermite (Brownian) and Laguerre 
(squared Bessel) cases A n = n. Hence, we will try to subordinate the Jacobi process by 
the mean of a random time-change in order to get a Mehler type formula. What is quite 
interesting is that subordinated Jacobi process semi-group, say qt(x,y), is the Laplace 
transform of P2/t{ x -iU)- Thus, we recover a suitable expression for pt(x,y) by inverting 
some Laplace transforms already computed by Biane, Pitman and Yor (see [I], |14|). 
This will allows us to derive a LDP for the maximum likelihood estimate (MLE) for p 
in the ultraspherical case, i. e. q = (/3 = a), a fact conjectured by Zani in her thesis. 
Then, using a skew product representation of the Jacobi process involving squared Bessel 
processes, we construct a family {vt}t of estimators for the index v of the squared Bessel 
process based on a Jacobi trajectory observed till time t. This satisfies a LDP with the 
same rate function derived for the MLE based on a squared Bessel trajectory. 

1.1. Inverse Gaussian subordinator. By an inverse Gaussian subordinator, we mean 
the process of the first hitting time of a Brownian motion with drift B^ := Bt+pt, p > 0, 
namely, 

Ip 5 = inf{s > 0; Bg = 5t}, 5 > 0. 
Using martingale methods, we can show that for each t > 0, u > 0, 

W e -» T t' S ) = e -t6(y/2u+»*-n) 

whence we recover the density below : 

= J^^exp 1 {S >0} 



^(P"' 13 (x)) n >o are normalized such that they form an orthogonal basis with respect to the probability 
measure W(y)dy which is not the same used in [T]. 
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1.2. The subordinated Jacobi Process. Let us consider a Jacobi process (Xt)t>o- 
Then, using (1), the semi-group of the subordinated Jacobi process (X T „,s)t>o is given 
by: 



Qt(x,y) = / p s (x,y)v t (s)ds 
Jo 

= W(y) ^(Rn)- 1 ( [ X e- x - s u t (s)d S ) P^(x)P^(y) 

n>0 ^0 J 



W(y)^2(R n r 1 E(e- x ^' S )P^(x)P^(y) 



n>0 

Writing A n = (n + 7) 2 — r y 2 where 7 = — — — , and substituting 5 = 1/V2, \i = V^T 
for a + (3 > — 1 in the expression of v t , one gets : 

so that 

q t (x,y) = W{y)Y J (Rn)- 1 e- nt P^{x)P^(y) 

n>0 

The last sum has been already computed ([T . p. 385) : 

1 - r ST^ \2) m+n \ 2 J m+n U V 



^0 ( x + r ) ™t>o ( Q + ^ + ^ mlnl 

1 — r „ .a a + 1 „ . 

(2) = (f^ * ( 2'~ ,a+ ;U,t,) 

where |r| < 1, a = a + /3 + 2, F4 is the Appell function and 

(1 — x)(l — y)r (1 + x)(l + y)r 

(l + r) 2 (1 + r) 2 ' 

Then, we use the integral representation of F4 (see |9j, p 51) to get: 

qt(x, y) = ^ nl\ a r s^e-'oFtia + 1; j S 2 )o*i(/3 + 1; ^s 2 )ds 
T{a) (1 + r) a J 44 

w w(y) i-r r^i.-y ^C*) n 2n 

r(o)(l + r)«y Jj(a + l^ + lV 



(2) 1 - r 



-y 7 r(2 " +a) - ^(z)a» 



r(a) (l + r)°^(a+l) n C9 + l) r 



where in (1), we used (see [12], p 214) 



Fi(c; w(l - r)/2) Fi(d; «;(! + r)/2) = ^ " ^ w n , a = c - 1, /3 = d-l, 



n>0 
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x ~\~ y (1 ~\~ xy)v 

in (2), we used Fubini Theorem, z = - — ; and A = — ; — — . Letting r = e~ t , then 



1 + xy 



_ W(y)e^ sinh(t/2) ^ 



(a) 



2/i 



2-i (cosh(t/2))« ^ (a + + 1), 



2(1 +r) 2 



(1 + xy) 
!cosh 2 (t/2) 



iy(y)tanh(t/2)e" 



Besides, 



2 a-l 



te7 t ,00 



E 

n>0 



(abn 



(a + !)„(/? + l) n 



1„ n 1 



(1 + xy) 



1 



cosh(i/2) 



2n+a-l 



2 t 

p s (x, y) s~ 3/2 e~ 7 s e~3i(is 



te7i ,00 



io 

Thus, noting that 7 = (a — l)/2, we get : 



P2/r{x,y)r 



P2/r{x,y)r 



-l/2 e -2 7 Vr e -Vr dr 



2^/2^.70 
27rW(y) tanh(t/2) 



-l/2 e -2 7 Vr e -V^ r 



E 

n>0 



(a) 



2/i, 



(a + !)„(/? + l) n 



P n ^(z) 



2 a-l 

1 + xy) 



t/2 

I \ 2n+a— 1 

cosh(t/2) y 



1.3. The Jacobi semi-group. The following results are due to Biane, Pitman and Yor 

(see [I], HU) : 



(3) 
(4) 



OO t 2 

e~^ s fc h ( s ) ds 







1 



00 2 



cosh(i/2) 
/tanh(i/2)\ h 



/i > 



ft > 



where (C/i) and (J/J are two families of Levy processes with respective density functions 
fc h and fx h for fixed h > 0. The densities of Ch and Ti are given by ([3|): 



2^l _1 r ^ Jr(2p+h){S) 



T(h) 



p>0 



M) = E e " V(fc+l)2si 



pi 



{s>0} 



k>0 



where r(c) = inf{r > 0; B r = c} is the Levy subordinator ( i. e, the first hitting time of 
a standard Brownian motion B) with corresponding density : 

(2p + h) ( (2p + h)^ 
fr(2 P +h){s) = 7= — = exp 



V27TS 3 



Thus : 

P2/r\X,y) - x ^ 



(a) 



2/i 



n>0 



(a +!)„(/? + 1), 



2s 



P n ^(z) 



L {s>0}- 



[l + xy) 



[fib */cw_i) (r 
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or equivalently (where B stands for the Beta function) : 



(a) 



2n 



2 a +f 3 Z_. (a + 1)n(/3+1)) 



il + xy) 



{fTl * fC2n+a + fi+l) (" 



1.4. The ultraspherical case. This case corresponds to a = (3 > — — and we will 
proceed slight differently. Indeed, a = 2a + 2 and 
1 — r 



(1 + r )2 Q +2 i?4 ( a + 1 ' a + 3 / 2 ' a + 1> a + 1; 



1 — r 



1 



(l+r) 2Q + 2 (l_ u _u)a+3/2--^ 4 ' | 

where the last equality follows from (see [5]) 



_ , 2a + 3 2a + 5 4ra , 
2 F i ( — ; — , — ; — ," + 1;^ — — ) 



F±(b,c,b,b;u,v) = (1 - u - v)- c 2 F 1 (c/2, (c + l)/2,b; 



(1 — u — v)< 



Hence, 



, . W(y)e^ lt sinh(t) ^,2a + 3 2a + 5 



(1 — « — u) 2 

(l-x 2 )(l-y 2 ) 



2 a+i/2 ( cos ht-xy)«+ 3 / 2 ' 



W(y)e 



2 a+l/2 



sin 



[(2a + 3)/4] n [(2a + 5)/4] n [(1 - x 2 )(l - y 2 )]' 



(cosht — xy)' 
2\m 



n>0 



(a + l)n 



(coshi - X y) 2n+a + 3 / 2 ' 



Besides, for h > 0, we may write : 



1 



cosh t — xy 



E 

fc>0 



(h) k (xy) h 



k\ (cosht) k+h 



since 



xy 



cosht 



< 1 Vx, y e] — 1, 1[, Vt > and where we used: 

7 -^- 7 - = V^r fe fc>0,|r|<l. 



fc>0 



Consequently, using Gauss duplication formula, 

T(u(n, k, a) + l){xy) k \{l - x 2 ){l - y 2 ) 



q t (x,y) = K a W(y)e^ t tanh{t) E 



n,k>0 



k\n\T(a + n + 1) 



1 



cosht 



v(n,k,a) 



where u(n,k,a) = 2n + k + a + l/2 and K a = T(a + l)/[2 a+1 / 2 r(a + 3/2)]. Thus, since 
7 = a + 1/2 when a = (3, one has : 



Ps{x,y)s 3/2 e rs e l*ds = -^ 



2vrr(a + 1) tanh(i) 



2 a r(a + 3/2) t 



W(y) 



E 

n,k>0 



T(v{n,k,a) + l){xyf 
k\n\T(a + n + 1) 



(l-x 2 )(l-y 2 ) 



2\l n 



1 



cosht 



u(n,k,a) 



6 
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or equivalently: 



00 _rl _£ 

Pi/2s(x, y)e 2s e 2 





ds ^T{a + 1) tanh(i) 
7s ~ 2 a r(a + 3/2) t 

T(u(n, k,a) + l)(xy) A; 



E 

n,fc>0 



(l-x 2 )(l 



fc!n!r(a + n + 1) 



1 



coshi 



u(n,k,a 



Using (j3J), ([3]), /c h et (we take i 2 /2 instead of t 2 /8), the density is written : 



Pi/2s{x,y) 



7rsT(a + 1) 



2°r(a + 3/2) 

ET(i/(n,k,a) + 1) fxy\ k 
k\n\T(a + n + 1) \~2~J 



Finally 

Pt(x,y) 



n,k>0 



W(y)e 
) 

) 



(l-x 2 )(l-y 2 ) 



/Ti * fc v r n k a s (s) 



E 

n,fc>0 



r(i/(n,ife,a) + l)(xy)* 
A;!n!r(a + n + l) 



(l-x 2 )(l-y 2 ) 



fax * fc v 



2. Application to statistics for diffusions processes 



2.1. Some properties of the Jacobi process. In the probability scope, we are used 
to define the Jacobi process on [—1, 1] as the unique strong solution of the SDE : 



dY, 



1 - Y t 2 dW t + (bY t + c)dt. 



It is straightforward that (Y t )t>Q = PQ/2)t>o where X is the Jacobi process already 
defined in section 1 with p = 26, q = 2c. In order to derive some facts, let us make the 
variable change y 1— > (y + l)/2, this gives up to a time change (t — ► At) : 



dJ t = 2v/J t (l - J t )dW t + [2(c - b) + AbJ t ] dt 



= 2y/j t (l - JtjdW t + [d-(d + d')J t ] dt 

where d = 2(c - b) = q - p = 2(/3 + 1) and d' = -2(c + b) = - (p + g) = 2(a + 1), which 
is the Jacobi process of parameters (d,d') already considered in |16j . Moreover, authors 
provide the following skew-product : let Z\,Z2 be two independent Bessel processes of 
dimensions d, d! and starting from z, z' respectively. Then : 



Z\(t) 



zf(t) + zi(t) 



(JA t ) 



t>0 ' 



t>o 



ds 



zf(s) + zi( s y 



Jn 



Using well known properties of squared Bessel processes (see |15j). one deduce that if 
d > 2 (/? > 0) and z > 0, then Jf > a. s. for all t > 0. Since 1 — J is still a Jacobi 
process of parameters (d',d), then, for d' > 2, (a > 0) and 2/ > 0, Jt < 1 a. s. for all 
t > 0. These results fit in the one dimensional case those established in [7] for the matrix 
Jacobi process (Theorem 3. 3. 2, p. 36). Besides, since is a reflecting boundary for 
Z\, Z2 when < d, d' < 2 (—1 < a, (3 < 0), then both and 1 are reflecting boundaries 
for J. 



LARGE DEVIATIONS FOR STATISTICS OF JACOBI PROCESS 



7 



2.2. LDP in the ultraspherical case. Let us consider the following SDE correspond- 
ing to the ultraspherical Jacobi process: 



(5) 



dY t = y/i - Y t 2 dW t + bY t dt 
Y = y E]-l,l[ 



Let Qy Q be the law of (Y t , t > 0) on the canonical filtered probability space (Q, (Tt),?) 
where Q. is the space of ] — 1, 1 [-valued functions. The parameter b is such that b < — 1 
(or a > 0), so that — 1 < Yt < — 1 for all t > 0. The maximum likelihood estimate of b 
based on the observation of a single trajectory (Y s , < s < t) under Qq is given by 

JO T^Yj Ub 

The main result of this section is the following theorem. 

Theorem 1. When b < —1, the family {bt}t satisfies a LDP with speed t and good rate 
function 

( l(x-b) 2 , 

t / \ I if X < Xq 

(7) J b (x) = I 4 x + 1 _ " 

[ x + 2 + ^/(b - x) 2 + 4(x + 1) i/ x > x > b 

where xq is the unique solution x < — 1 of the equation 

(b - x) 2 = 4x(x + 1) . 

Proof of Theorem [J^ 

We follow the scheme of Theorem 3.1 in [19J. Let us denote by : 

rt Y s f f Y 2 
Six = / — rr dY s — XI - — k ds 

' Jo 1~Y S 2 Jo l-Y 2 

so that for x > b (resp. x < b), P(b t > x) = P{S t , x > 0) (resp. P(b t < x) = P{S t , x < 0)). 
Therefore, to derive a large deviation principle on {bt}, we seek a LDP result for St iX /t 
at 0. Let us compute the normalized cumulant generating function At x of St iX : 

(8) A f = i log £(e<^) 
From Girsanov formula, the generalized densities are given by 

dQ 

From Ito formula 



1 rt y i rt i I y2 

Let us denote by 

V 1 = {<j> : (b + l) 2 + 2(j)(x + 1) > 0} . 
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Set b(</>, x) = -1 - ^/(b + 1) 2 + 2(f){x + 1) for all <\> G Di. Then : 

A t (<f>,x) = hogE b{(j)jX) (eM{<t> + & " b(<f>,x))[F(Y t ) - F(y ) - t/2]}) 
When starting from yo = 0, the semi- group is deduced from that of X : 

n>0 

where p = -2(a + 1) = 2b < -2 and 7 = -(p+ l)/2 = a + 1/2. Denote by 

P = {^GPi: GO) = 6 + 60, x) + < 0} . 
For any <fi £ T>, the expectation above is finite and a simple computation gives : 

MM = - 4> + b - b ^ x) + %iW(i " i?r<"^*»/ 2 ) 



= A(0, x) + - log ^ 

where 

a W x) = V r(2n-6Q,x) + l/2) / _ + 6 + 1\ 2t/ 1 

and B stands for the Beta function. Moreover, by dominated convergence theorem 
lim E 6W>X) ((1 - y t 2)-(*+6-6(^))/2) = /' (1 - 2/ 2 ) -[^ +6+b (^)]/2-i dy < ^ 

t— >oc ' 

for <fi £ T>. Hence At — > A as i — > 00. The following lemma details the domain P of At: 
Lemma 1. Denote by 

_ (6+1) 2 

00 



2(x + l) ' 

i/ ffx < (b 2 + 3)/2(fe- 1) : iften D = (-00, <fo). 

raj // (b 2 + 3)/2(6 - 1) < x < -I: then V = (-00, (pi) where (pi is solution of G {(f) = 0. 

Hi) If x > — 1 : then V = (<po,4>i)- 

In case i) of Lemma above, A is steep. It achieves its unique minimum in <p m solution 

of 

dA/d<f>(<i>,x) =0, 
i.e. b(4>, x) = x. It is easy to see that 

x + 1 (b + l) 2 ^ 

-^TT) <<t>0 - 

Hence, Gartner-Ellis Theorem gives for x < b < (b 2 + 3)/2(& - 1), 

lim y log P(b t <x) = lim - log P(S t , x <0)= inf A(0, x) = A(^ m , x) = -\ (x ~ *f . 
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If b < x < (b 2 + 3)/2(& - 1), notice that <p m > and 



i) 



1 1 1 (x — bY z 

lim - log P{b t >x) = lim - log > 0) = inf A(<£,x) = A(0 m ,x) - " " 



t->oo i t-+oo i ' 0e (o,(^ ] 4 x + 1 

In cases ii) and hi) of Lemma [H A is not steep. Nevertheless, if the infimum of A 

o 

is reached in D, we can follow the scheme of Gartner-Ellis theorem for the change of 
probability in the infimum bound. This infimum is reached if and only if 

(9) 8A/d(j)((j)i,x)>0. 

This above condition gives the following cases: denote by xo the unique solution x < 
— 1 of g(x) := 4x(x + 1) — (b — x) 2 = 0. Since g is decreasing on ] — oo, — 1] and 
g(b 2 + 3/(2(6 - 1)) = (3/4)(6 + l) 2 > = g(x ), then x > (b 2 + 3)/ [2(6 - 1)]. 

• if (b 2 + 3)/2(6 — 1) < x < xq < — 1, the derivative dA/d<j)((f>i,x) > 0, A achieves its 
minimum on <j) m and 

1 (x - b) 2 

lim -logP(6 t > x) = A(cpm,x) = - \ ' . 
t->co t 4(cc + 1) 

• if xq < x < — 1 or x > —1, then 8A/d(f)(<f)\,x) < 0. We apply Theorem [2] of the 
appendix, which is due to Zani [19]. Let us verify that the assumptions are satisfied. 
Indeed, the only singularity <j>\ of Rt comes from B (n — [eft + b + &(</>, x)]/2, 1/2) when 
n = 0, and more precisely, from r (—[</> + b + b(<j), x)]/2) . We can write 

(10) A t (<j),x) = A(4>,x) + -logT + - log- 



where 

(U) *(* I > = r(-[* + » + 6(*,x)]/2) 

Now 

/ </> + 6 + 6(</>,x) 1 

£> n , - 

V 2 2 

Vn > 0, — > 



r 



+ 6 + 6(0,x) 



2 

is analytic on some neighbourhood of <j>\. Besides, <j>\ is a pole of order one, i.e. 

r b + (f> + b(cf>,x) 

hm = c > , 

4>^4>i ,4><<t>\ 4> — <pi 

and since lim p ^ + p£(p) = 1, we can write 

1 / ( f ) + b + b(4>,x) \ log(^-0) h{<p) 

t og V 2 / 1 1 

The function h is analytic on D and can be extended to an analytic function on \<f)\ — 
£> 4>i +£[ f° r some positive ^. Finally, we focus on Rt{4>, x) / \/i and show that it converges 



10 



LARGE DEVIATIONS FOR STATISTICS OF JACOBI PROCESS 



uniformly as t — > oo. To proceed, we shall prove that this ratio is bounded from above 
and below away from 0. Setting A n (t) := e 7 t ^ 2 fT 1 */c 2 „+ 7 one nas • 



A n (t) < e 



y/t " y/t 
P ft/2 



1/t 1 

exp — - 



(2„ + 2 t + 7 )3 +j2(i l )2( l s) 

S Zi Tj 



E u ^ f 

k,l>0 J ° 

/oo i j 

exp-- (2n + 2k + 1 fs + K 2 {l + -f{ 
« ... L 



l N o,s — t s 



ds 

13/2 



/OO -1 
exp-i 



is J \A 



( 2n + 2k + 7 ) 2 (s - 1) + vr 2 (/ + i) 2 ( i-^) 



2' v is 'J Vte 



cis 



k,l>0 



2/r 2 / 1 



(2n + 2/c + 7) 2 s + vr 2 / 2 ( 



t(t + S ) J \ y/i(t + 8) 



with 



T{2n + k + 7)2 2 " +7 (2ra + 2k + 7) 
/c!r(2n + 7 ) 



Let G(x) = Y,i& e x = 1 + 2 ^>i e _7r ' * denote the Jacobi Theta function. Then 



A n (t) 

where 



< e 



-2n" 



_^ /-oo 

/ exp- 
fc>0 Jo 



(2n + 2/c + 7)V 



e 



7TS 



(is 



1 2 r 00 

C(n, i) = — = V U Kn e~ 2k2 / exp ■ 

2 ^ fc,i>0 7o 



2i(i + s)y 7i(iT7 

(2n + 2A; + 7) 2 sl cis 



+ C(n,t) 



v^Ts 



Recall that 6(x) = (l/y/x)Q(l/x), which yields : 



where 



< e" 2 " 2 E U Kn e~ 2k2 / exp ■ 



(2n + 2/c + 7) 2 s 1 



e 



2i(i + s)^ cis C(n) 



C(n) = e- 2 " 2 ^ C/ fc , n e- 2fe2 / 



exp 



7rs / y^s 2^/i 
(2n + 2/c + 7) 2 sl ds 



Since e _/ z < e~ lz , then B(z) < 3 for z > 1. Hence, as 2i/7r < 2i(i + s)/(ns), then for 
i large enough : 



A n {t) 

Vt 



„ /-oo 

< 3e- 2 " 2 ^ [/ M e" 2fc2 / exp 



(2n + 2k + 1 y S ds_ 
2 
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This gives a lower bound for R±/\Jt. Besides, 



R t (<j),x) 



> 



v^r(l/2-6(0,x)) 



A (t) 



yft r(-6(0,x))r{[l-(0 + 6 + 6(0,a:)]/2} ft 



exp-- 

/,-./ ii 



2' t(t + s) 



ds 



s/t{t + s) 



where V k (t) := C/ fc ,oe" 2fc(fc+7) *. One may choose t large enough independent of k such 
that Vfc(t) > Vfc + i(i) for all k > 0. In fact, such t satisfies : 

e 2( 2fe+7+ i) t > g2t > gup £W = gup .(* + 7) (2* + 7 + 2) 



fc>0 



t4,o fc>o (A; + l)(2fe + 7) 



Then 



-|>C(6,0,x)[F o (t)-^(t)]^ / exp-- 

/■OO -1 

>C(6,0,x)[ 7 2^-y 1 (t)]^ / exp-- 

Z>0 ^° 

C(6, <p, x 



7 2 » + » 2 (I+l) 2 ( " 



ds 



t(t + a) 



ds 



-7 2 s/2@ 



ITS 



ds 



2t{t + S )J y/iit + sj 



C(t) 



where 



C(t) 



1 



for t large enough. Following the same scheme as for the upper bound, one gets 



2V~tJ 
Follow 

Vt > 2 



°° -y 2 s/2 



< C 



72 7 < \ ~ 



e _ 7 2 s/2@ (2t(t + s)\ ds 



ITS J y/s 



-= - C(t) 



> 



C(b,<j>,x) 



7 2T W-- c 



□ 



As a result, 



lim -logP(6 t > x) = A(0i,x) = -(x + 2 + */(& - x) 2 + 4(sc + 1)), 

which ends the proof of Theorem [H ■ 

2.3. Jacobi-squared Bessel processes duality. By Ito's formula and Levy criterion, 
one claims that (Y t 2 ) t >o is a Jacobi process of parameters d = 1, d! = —2b > 2. Indeed : 



dZ t := d{Y t z ) = 2Y t dY t + (Y) t = 2Y t Jl - Y t 2 dW t + [{2b - l)Y t z + l]dt 



2yJZ t {l - Ztjsgn(Y t )dW t + [(26 - l)Z t + l}dt 



2 V / Z t (l - Z t )dB t + [(26 - \)Z t + l]dt 
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Using the skew product previously stated, there exists R, a squared Bessel process of 
dimension d! = 2{u + 1) = —2b and starting from r so that : 

u n iog(i-y t 2 )+* 



2 /o w?* 



is another estimator of v based on a Jacobi trajectory observed till time t. Set t = log-u, 
then 

_ log[u(l - Y^J] _ log[«(l - Y£ S J] 



1A 



log" r log« y s 2 , y2 

V logs J 

and {^i 1 ogn }u satisfies a LDP with speed logu and rate function J_(j, +1 )(— (x + 1)). 

When starting at Rq = 1, the MLE of v based on a Bessel trajectory is given by (cf 
M- P- 132) : 

*1 = So ^YT ~ 2 Jo it = log TO 
z Jo X s z Jo X s 

with associated rate function : 



i Ax) 



(x-v) 2 :f _ > - ._ -(^+2)+2yV2 +;y+ i 



1 — x + \J (y — x) 2 — 4x if x < a?i 
A glance at both rate functions gives I u (x) = J_(V +1 )(— (x + 1)) and xq = —{x\ + 1). 

3. Appendix 

Let {Yt}t>o be a family of real random variables defined on (O, J-, P), and denote by 
/it the distribution of Y t . Suppose — oo < rat '■= EY t < 0. We look for large deviations 
bounds for P(Y t > y). Let A t be the n.c.g.f. of Y t : 

A t (0) = ilog E(exp{<f>tY t }) , 

and denote by D t the domain of A*. We assume that there exists < <f>\ < oo such that 
for any t 

sup{(j) : <p G D t } = 4>i 
and [0, 4>i ) C Df. We assume also that for <p G D 

Assumption 1. 

(12) A t (0) = A(<£) - f log(0i - 0) + ^ 

where 

• a > 

• A is analytic on (0, (f>\), convex, with finite limits at endpoints, such that A'(0) < 0, 
A'(0i) < oo, and A"Oi) > 0. 

• Rt is analytic on (0, cpi) and admits an analytic extension on a strip Z}^ = {<pi — /3,cpi + 
(3) x (—7,7), where (3 and 7 are independent of t. 

• Rt((p) converges as t — > 00 to some i?(^>) uniformly on any compact of DZ. 
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Theorem 2. Under 

For any A'(0) < y < A' (fa), 

(13) lim - log P(Y t > y) = - sup {y</> - A(<j>)} . 
For any y > A' (fa), 

(14) lim \\ogP(Y t >y) = -yfa + A(fa). 
The rate function is continuously differentiable with a linear part. 
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